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ABSTRACT 

We develop a semi-analytic model for planet formation during the pre-main 
sequence contraction phase of a low mass star. During this evolution, the stel- 
lar magnetosphere maintains a fixed ratio between the inner disk radius and 
the stellar radius. As the star contracts at constant effective temperature, the 
'snow line', which separates regions of rocky planet formation from regions of 
icy planet formation, moves inward. This process enables rapid formation of icy 
protoplanets that collide and merge into super-Earths before the star reaches the 
main sequence. The masses and orbits of these super-Earths are consistent with 
super-Earths detected in recent microlensing experiments. 



Subject headings: planetary systems: formation — planetary systems: proto- 
planetary disks — stars: evolution — stars: formation 
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1. Introduction 

The recent discoveries of super-Earths around low mass stars challenge our understand- 
ing of planet formation. From separate microlensing events, Beaulieu et al. (2006) and Gould 
et al. (2006) provide strong evidence that planets with masses M ~ 5-15 M e are common 
around M dwarf stars 1 . With orbital semi-major axes a ~ 2.5-3 AU, these planets are 
probably ice giants roughly similar in structure to Uranus and Neptune in the Solar System. 

Boss (2006) proposes that these planets form in two stages. After a disk instability 
produces a gas giant, photoevaporation of the gas giant atmosphere leaves behind an icy 
core with M ~ 10-20 M e . This mechanism requires a massive disk to initiate the instability 
and a nearby O-type star to photoevaporate the gas giant atmosphere. Boss notes that this 
process should yield (i) super-Earths around M dwarfs formed in rich star clusters and (ii) 
gas giants around M dwarfs formed in low mass stellar associations. 

Beaulieu et al. suggest that super-Earths favor coagulation models, where collisions of 
1-10 km objects eventually produce icy planets with M ~ 10 M ffi at 1-10 AU. Although 
numerical calculations appear to preclude gas giants at 1-10 AU around M dwarfs (Laughlin 
et al. 2004), there has been no demonstration that coagulation produces icy planets on 
reasonable timescales in a disk around an M dwarf. 

Here, we develop a semi-analytic coagulation model, and show that contraction of the 
central star along a pre-main sequence (PMS) Hayashi track sets the initial conditions for 
planet formation around low mass stars. Our results indicate that icy protoplanets with M ~ 
0.1-1 M e form in ~ 0.1-1 Myr at 1-4 AU. Over 50-500 Myr, collisions between protoplanets 
produce super-Earths with masses similar to those detected in microlensing surveys. 

We start with the motivation for our study in §2, discuss the coagulation model of planet 
formation and the moving snow line in §3, develop the disk evolution model in §4, and apply 
the model to super-Earth formation in §5. We end with a brief summary in §6. 

2. Motivation: planet formation in the disk of a low-mass star 

To motivate our study, we contrast planet formation around low mass stars and solar- 
type stars. For solar-type stars approaching the main sequence, the luminosity is roughly 



1 These microlensing observations directly yield the ratio Mp/M±, where Mp is the mass of the planet 
and M* is the mass of the star. Adopting the most likely parent star - a K-dwarf or an M-dwarf - yields 
the most likely Mp. 
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constant on typical planet formation timescales of 10-100 Myr. Thus, the conditions where 
planets form change little with time. For stars with masses < 0.5 M , however, the lumi- 
nosity fades by a factor of 10-100 on the Hayashi track. Because the inner disk radius is 
'locked' at a fixed distance relative to the radius of the central star, the inner disk contracts 
as the star contracts. During this evolution, the 'snow line' - the point that separates the 
inner region of rocky planet formation from the outer region of icy planet formation - also 
moves inward. 

Coupled with the evolution of the inner disk, the moving snow line produces a dramatic 
variation in the surface density at fixed distances from the central star. This behavior 
enables the rapid formation of icy protoplanets. As the low mass star approaches the main 
sequence, these protoplanets collide and merge into super-Earths with properties similar to 
those detected in recent microlensing experiments. 



In coagulation models, planets grow from repeated collisions and mergers of small objects 
in a circumstellar disk (Safronov 1969). When 1-10 km 'planetesimals' form and start to grow 
(Weidenschilling 1980; Dullemond & Dominik 2005), dynamical friction damps the orbital 
eccentricities of the largest objects. Damping yields large gravitational cross-sections and 
leads to 'runaway growth,' where the largest objects grow fastest and run away from more 
slowly growing smaller objects (Wetherill & Stewart 1989; Kokubo & Ida 1996). Throughout 
the runaway, the largest protoplanets stir up the leftover planetesimals. Eventually, the 
leftovers have orbital velocity dispersions comparable to the escape velocities of the largest 
protoplanets. Because gravitational cross-sections fall as velocity dispersions rise, runaway 
growth ends. The ensemble of planetesimals and protoplanets then enters 'oligarchic' growth, 
where the largest objects - oligarchs - accrete at rates roughly independent of their size 
(Kokubo & Ida 1998). 

During oligarchic growth, protoplanets become isolated from their surroundings. If an 
oligarch accretes all of the mass in an annulus with width 2BRh, where Rh = a(M/3M i ) 1 / 3 
is the Hill radius, its isolation mass is 



where a is the orbital semi-major axis and a is the mass surface density of solid material in 
the disk (e.g. Lissauer 1993; Kokubo & Ida 2000). If p is the mass density of a solid object 
the timescale to reach isolation is (Goldreich et al. 2004) 



3. Coagulation and the Moving Snow Line 



M iso « 4iraBR H a oc (Ba) 3/2 a 3 M- 1/2 , 




t lso oc pVV/V-1/2 



(2) 
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In the Solar System, oligarchic growth has two regimes. For a < 3AU, planetesimals are 
rocky because volatile materials remain in the gas. In the 'Minimum Mass Solar Nebula' 
(MMSN, Weidenschilling 1977; Hayashi 1981) with B = 2.5 - 5, p ~ 3 gem" 3 , and 
er ~ 8gcm~ 2 at 1AU, M iso ~ 0.05-0.1 M e . Once oligarchs contain ~50% of the total 
mass in solids, their mutual dynamical interactions lead to 'chaotic' growth (Goldreich et 
al. 2004; Kenyon & Bromley 2006), where collisions between oligarchs eventually produce 
Earth-mass planets. Numerical simulations suggest that ~10-20 oligarchs collide to form 
a typical Earth-mass planet in ~10-100 Myr at 1AU around a solar-type star (Chambers 
2001; Raymond et al. 2004; Kenyon k Bromley 2006). 

Outside the 'snow line', ice condensation enhances a and promotes the formation of 
larger oligarchs. For p ~ 1.5 g cur 3 and a ~ 3-6 gem -2 at 5AU, isolated oligarchs with 
M iso ~ 5 M e form on timescales t iso ~ 1-3 Myr. These icy oligarchs accrete gas directly from 
the nebula and grow into gas giant planets in several Myr (Pollack et al. 1996), comparable 
to the lifetime of the gaseous disk (e.g. Hollenbach et al. 2000; Haisch et al. 2001; Young et 
al. 2004; Calvet et al. 2005). 

For solar-type stars, planet formation is fairly independent of stellar evolution. Through- 
out most of the PMS phase, the solar luminosity is roughly constant. Thus, the position of 
the snow line - a snow ~ (L^/T^^) 1 / 2 , where L+ is the stellar luminosity and T snow is the 
temperature where water and other volatile materials condense out of the gas - is roughly 
stationary in time. Because the ~0. 1-1 Myr formation time for planetesimals and oligarchs 
is short compared to the ~10Myr PMS lifetime, the separation between icy and rocky 
(proto)planets remains fairly distinct, evident in the composition of different populations 
in the asteroid belt (Abe et al. 2000; Rivkin et al. 2002). Although there is some mixing 
between water-rich and water-poor regions (Raymond et al. 2004), most of chaotic growth 
occurs when the Sun lies close to the main sequence at nearly constant L+. 

In contrast with solar-type stars, stellar evolution is a crucial feature that defines the na- 
ture of newly-formed planets around low mass stars. Because the timescale for planetesimal 
and oligarch formation is short compared to the 0.1-1 Gyr PMS contraction time (D'Antona 
& Mazzitelli 1994; Baraffe et al. 1998; Siess et al. 2000), the timing of planetesimal forma- 
tion sets the nature of icy/rocky planets with distance from a low mass star. On its Hayashi 
track, the luminosity of a 0.25 M & star fades by a factor of several hundred at roughly con- 
stant effective temperature. During this period, a snow moves inward by a factor of ~15-20. 
Just outside the moving snow line, ice condensation increases a (M iso ) by a factor of ~ 4 
(8) (Hayashi 1981); t iso decreases by a factor of 3. This moving snow line enables rapid 
formation of icy oligarchs that can collide and merge into super-Earths. 
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4. Evolution of a disk around a contracting star 



Disk evolution is also an important feature of planet formation around low mass stars. 
In the standard MMSN model, a is fixed in time and scales with the stellar radius on 
the main sequence (e.g. Hayashi 1981). However, when PMS stars actively accrete from a 
circumstellar disk, magnetic interactions between the star and the disk appear to 'lock' the 
inner disk radius Ri n at a fixed distance relative to the stellar radius, £ = Ri n / R+ ~ 3, 
at several Myr (e.g. Eisner et al. 2005). Although the duration of this phase is not well- 
constrained, the observed change in £ for disks around solar-type stars is a factor of ~ 2-3 
(Eisner et al. 2005). If disks around low mass stars remain locked for the entire PMS phase, 
the maximum decrease in the inner disk radius is a factor of ~15-20. This change is much 
larger than the observed variation of £; thus we assume £ = constant. To conserve mass and 
angular momentum, a and the outer disk radius must evolve, which impacts M iso and the 
formation timescales for oligarchs and planets. 

To construct a model for disk evolution, we adopt 



where R* is in units of solar radii, <r = 8 gem -2 , and a AU is the radial distance from the 
star in AU. Setting the scale factor (3 ~ 3 yields the usual cr(MMSN) for a 1 M star at 
lMyr, when a large fraction of the solid mass in the terrestrial zone of the Solar System 
is in large bodies. Consistent with observations (Natta et al. 2000; Scholz et al. 2006), we 
scale a and the disk mass linearly with the stellar mass. For a lMyr old 0.25 M & star, this 
scaled MMSN has (3 = 2 and Mdisk = 0.026M* integrated from 3.R* to 50 AU for a gas/solids 
ratio of 100. To provide a smooth transition from f ice = 1 for a < a snow to fi ce = 4 for 
a ^ a snow (Hayashi 1981), we include a parameter f ice = 1 + (A ice — 1)/(1 + e x ) where 
A ice = 4, x = (a snow — a) / AT snow (a) and AT snow (a) is the radial distance equivalent to a 5 K 
temperature change. 

To derive a snow , we adopt the temperature profile of a flat circumstellar disk, T oc 
(i?*/a) 3//4 (Kenyon & Hartmann 1987). We scale this relation to place the snow line at 
2.7AU at lMyr for a 1 M mass star, as inferred from analyses of water-rich objects in 
the outer asteroid belt (Abe et al. 2000; Rivkin et al. 2002). To evaluate L±(t), R+(t), and 
T+(t), we use PMS evolutionary tracks from Siess et al. (2000); other tracks (D'Antona & 
Mazzitelli 1994; Baraffe et al. 1998) yield similar results. 

With these ingredients, we derive the evolution of a, M iso , and t iso as the star contracts 
to the main sequence. This evolution has two main features. Initially, the snow line is at a 
large distance, a snow ~ 5AU, from the luminous PMS star. Inside 1-2 AU, rocky oligarchs 




(3) 
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form and reach M iso before the star contracts significantly. Outside 1-2 AU, t iso is long 
(eq. 2) compared to the initial contraction time. As the star contracts, ices condense out of 
the nebula and the snow line moves inward. For a < 1-2 AU, this material coats the growing 
oligarchs, leftover planetesimals, and the surrounding debris with an icy veneer that may 
extend the oligarchic growth phase and produce more massive oligarchs. For a > 1-2 AU, 
ice condensation reduces ti SO by a factor ~ 3 (eq. 2), which enables the rapid formation of 
icy oligarchs well before the central star reaches the main sequence. 



5. Super-Earth Formation 

To explore the consequences of this picture, we consider a disk with /3 ~ 1 (M^^/M* = 
0.065), which lies at the upper end of the range inferred from observations 2 (Osterloh & 
Beckwith 1995; Nuernberger et al. 1997, 1998; Natta et al. 2000; Scholz et al. 2006). Figure 
1 shows the a evolution for this system at several distances from a 0.25 M Q star. For disks 
with other masses, a oc M disk , M iso oc M^/ S 2 fc , and t iso oc M^J k 2 . Aside from the long-term 
decline in a (a) from PMS evolution, the a evolution shows clear increases when the snow line 
crosses specific points in space and ices condense out of the gas. At these times, a remains 
at a relatively constant plateau value for ~l-3Myr before declining monotonically as the 
central star approaches the main sequence. 

During the plateau phases, the timescale for oligarchs to reach M iso (< 1 Myr; eq. 2) 
is shorter than the stellar contraction time. Because the isolation masses at fixed distances 
decrease as a decreases (eq. 1), these early times provide the best opportunity for coagulation 
to form large protoplanets. 

Figure 2 shows the time evolution of M iso . Interior to a snow (a < 1-2 AU), rocky oligarchs 
with Mi SO ~ 0.1 Mq form in ~10 5 yr. As the star contracts, ice condensation enables the 
formation of larger oligarchs with M iso ~ 0.2 M e in ~1 Myr. At a ~ 2-3 AU, ice condensation 
during runaway growth promotes the formation of oligarchs with M iso ~ 0.5 M e in ~ 10 5 yr. 

This analytic prescription for protoplanet growth suggests that oligarchs with M iso ~ 
0.1-0.5 M e can form at ~1-3AU in < 1 Myr. If oligarchs contain roughly half the mass in 
solid material at the onset of chaotic growth, our model disk with ~ 5-10 M e at 1-4 AU will 
have ~10-100 oligarchs. The model predicts ~10 oligarchs at 2-4 AU. Thus, the building 
blocks for observable super-Earths can form on timescales much shorter than disk lifetimes 



2 In their coagulation model for Neptune, Goldreich et al. (2004) also consider a disk with Mdisk ~ 3- 
GMmmsn- 
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derived from measurements of dust emission from low mass PMS stars (Song et al. 2002; 
Weinberger et al. 2004; Liu et al. 2004; Plavchan et al. 2005). 

To consider whether oligarchs can merge into super-Earths on reasonable timescales, 
we follow Goldreich et al. (2004) and introduce a parameter 1Z = v esc /Qa, where v esc is the 
escape velocity and ila is the orbital velocity of a growing planet. When K« 1, colliding 
oligarchs merge; when 1Z ^> 1, collisions often eject an oligarch from the planetary system. 
At 1-5 AU around a 0.25 M star, this merger condition (1Z < 1) allows the formation of 
~5M ffi planets at 1-2 AU, ~ 3M e planets at 2-3 AU, and 1-2 M© planets at 3-4 AU. In 
this analytic model, the timescale to produce planets is ~ 200 Myr (1 Gyr) at 1 AU (5 AU). 

To derive another estimate for the masses and formation timescales, we consider the 
results of complete numerical simulations of planet formation from an initial ensemble of 
oligarchs (Chambers 2001; Raymond et al. 2004) or planetesimals (Kenyon & Bromley 2006). 
In the solar terrestrial zone, collisions and mergers of 10-20 oligarchs with masses ~ Mi so 
yield 2-5 planets with masses comparable to the mass of the Earth on timescales of 10- 
100 Myr. Although the final orbital parameters depend on the late-time evolution of the 
planetesimals and the gaseous disk, the typical masses and collision histories of Earth-mass 
planets are similar in all calculations and agree fairly well with analytic estimates. Adapting 
this collisional history to a planetesimal disk around a 0.25 M Q star, mergers of ~ 10 oligarchs 
should yield planets with masses ~l-2M e at 1 AU and ~ 3-5 M e at 2.5 AU. 

Combining the analytic and scaled numerical results, the timescale for oligarchs to merge 
into planets is roughly 



where P is the orbital period. Thus, the expected merger timescale for oligarchs at 1-3 AU 
around a 0.25 M star is ~2-5 times longer than for the terrestrial zone around a solar- 
type star. This timescale is comparable to the duration of the PMS contraction phase and is 
much shorter than the expected stellar lifetime. Thus, coagulation can produce super-Earths 
around low mass stars on timescales of ~ 50-500 Myr. 

In constructing our model, we adopted a standard surface density law, a oc a~ 3 / 2 , 
and ignored details of the disk structure (e.g. Davis 2005; Lecar et al. 2006) and physical 
mechanisms for ice condensation (e.g. Podolak & Zucker 2004). Although details of the disk 
structure and ice condensation mechanisms can affect the position of the snow line, our main 
conclusions that (i) a snow moves considerably during the PMS contraction of a low mass star, 
and (ii) ice condensation during the PMS contraction phase produces massive oligarchs in 
~ 0.1-1 Myr and super-Earths in ~ 100 Myr, are generally independent of these details. 
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The main uncertainties in our picture are the probability of the large initial disk mass and 
the details of the final accretion stage when 1-2 M e planets evolve into 3-5 M e planets. 
Observations of larger samples can yield better estimates for the range of initial disk masses 
for low mass stars and for the M^fe-M* relation. Detailed numerical simulations can provide 
better estimates of the masses and formation timescales for super-Earths. 

6. Model summary and predictions 

We have developed an analytic prescription for planet formation by coagulation around 
low mass stars. The model has two distinctive features that enable formation of super-Earths 
during the PMS contraction phase. 

• We set the inner disk radius at a fixed distance relative to the radius of the central star, 
£ = R in /R±. Thus, R in changes as the star contracts to the main-sequence, leading to 
significant evolution in a (a). 

• During PMS contraction, a snow moves inwards by a factor of ~ 15-20, producing large 
enhancements in a (a) as ices condense out of the nebula. Ice condensation is the key 
mechanism that allows coagulation to produce super-Earths around low mass stars. 
This process results in new outcomes for planet formation, including planets with 
dense, rocky cores and thick, icy surfaces. 

We applied this model to super-Earth formation around a 0.25 M & star. At 1-5 AU, 
isolated oligarchs can grow to masses ~ 0.1-1 M e in ~ 0.1-1 Myr, short compared to the 
~ 100 Myr PMS contraction time. These oligarchs merge into super-Earths with masses ~ 2- 
5 M e as the star contracts to the main sequence. Thus, coagulation can produce planetary 
systems similar to those detected in recent microlensing events. 

Aside from our success in producing icy super-Earths at 1-3 AU around low mass stars, 
the model makes clear predictions for rocky planet formation close to a low mass star. At 
a < 1 AU, the isolation masses estimated for rocky planets are M iso < 0.01 M e . If ~10 
isolated objects merge into a rocky planet, we predict many low mass planets with masses 
~0.1M e at distances of 0.05-0.5 AU around low mass M dwarfs. We plan to explore the 
consequences of our model in future papers. 

We acknowledge support from the NASA Astrophysics Theory Program through grant 
NAG5-13278 (SK, BB) and an Australian Postgraduate Award (GK). We thank T. Currie, 
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Fig. 1. — Surface density evolution at fixed radii around a 0.25 M Q star with M disk jM+ = 
0.065. As the snow line moves inwards, ice condensation increases a, which leads to faster 
formation of more massive oligarchs. 



-14- 




Fig. 2. — Evolution of M iso at fixed radii using the a evolution of Figure 1. Ice condensation 
leads to more massive oligarchs at 2-8 AU in 1 Myr. 



